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ABSTRACT 

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the 
structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We 
included relevant processes to describe mass transfer between the different fluids, ra- 
diative cooling by emission lines and grain charging and studied the effect of single- 
sized and multiple sized grains on the shock structure. Our models are the first of 
oblique fast-mode molecular shocks in which such a rigorous treatment of the dust 
grain dynamics has been combined with a self-consistent calculation of the thermal 
and ionisation structures including appropriate microphysics. At low densities the 
grains do not play any significant role in the shock dynamics. At high densities, the 
ionisation fraction is sufficiently low that dust grains are important charge and current 
carriers and, thus, determine the shock structure. We find that the magnetic field in 
the shock front has a significant rotation out of the initial upstream plane. This is 
most pronounced for single-sized grains and small angles of the shock normal with 
the magnetic field. Our results are similar to previous studies of steady C-type shocks 
showing that our method is efficient, rigorous and robust. Unlike the method employed 
in the previous most detailed treatment of dust in steady oblique fast-mode shocks, 
ours allows a reliable calculation even when chemical or other conditions deviate from 
local statistical equilibrium. We are also able to model transient phenomena. 
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1 INTRODUCTION 

Molecular outflows associated with proto-stellar objects (e.g. 
iBallv fc Ladalll983D drive shocks into dark molecular clouds 
where the fractional ionisation \ is low (x < 10 -6 ). The low 
values of \ cause the gas and the magnetic field to be weakly 
coupled. Charged particles are then pushed through the neu- 
tral gas by Lorentz forc es, giving rise to ambipolar diffusion 
|Mestel fc Spitzerl 19561 ) . Collisions between the charged and 
neutral particles accelerate, compress and heat the neutral 
gas ahead of the shock front and a magnetic precursor forms. 
If the shock speed is low enough for the gas to cool efficiently , 
the flow becomes continuous in all fluids fe.g. lMullan| | l971f) . 
Suc h shocks are referred to as C-type sho cks (|Drainell 19801 ). 

iDraine. Roberge. fc Dalgarnol (Il983l . hereafter DRD) 
calculated the structures of steady perpendicular C-type 
shocks in diffuse and dense molecular media using a steady- 
state multifluid magnetohydrodynamic (MHD) description. 
They included a detailed treatment for energy and mo- 
mentum transfer between different particles, oxygen-based 
chemistry and radiative cooling. They also assumed that the 
ion and electron densities exceed the grain charge density. 

* E-mail: physvl@leeds.ac.uk 



However, for a fractional ionisation below 10 -7 , grains play 
an important role in governing the electromagnetic forces 
l|Drainel 1 1980T ) . In the densest regions of molecular clouds 
where ng ^ 10 6 cm -3 , the fractional ionisation constraint 
is usually fulfilled. 

IPilipp. Hartquist. fc Havnesl l|l990t ) included the effect 
of dust grains on the structure of perpendicular steady C- 
type shocks more rigorously than previous authors (e.g. 
DRD). Their four- fluid models predict different shock widths 
than obtained in similar mo dels of DRD. This model w as ex- 
tended to oblique shocks bv lPilipp fc Hartquisd (|l994l ) who 
found that in such shocks the tranverse component of the 
magnetic field rotates about the propagation direction of the 
shock as long as there is a non-negligible Hall conductivity. 
This happens for a wide range of conditions, but Pilipp & 
Hartquist were unable to find solutions for fast-mode shocks. 

IWardlei (1 19981 ') pointed out that, to find a stable shock 
structure for steady oblique fast-mode shocks, one needs to 
integrate the steady state equations backwards in time from 
the downstream boundary. However, such a method requires 
an exact thermal, chemical and ionisation equilibrium. Such 
an equilibrium does not obtain in most shocks of interest. 
For instance, the abundance of water, an important coolant, 
is out of equilibrium for about 10 5 years or more from the 
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time the gas cools below several hundred degree. This high- 
lights the necessity for time-depe ndent simulations t o cal - 
culate the C-type shock structure. ICiolek fc Robergd (|2002h 
were the first to include a fluid description of grain dynam- 
ics in multifluid MHD simulations. However, their study is 
restricted to perpendicula r shocks. The multifluid approach 
developed by Valid (|2003l ) overcomes this limitation. 

In this paper we study the structures of steady fast- 
mode C-type shocks using time-dependent multifluid MHD 
simulations. In Sect. [2] we describe the numerical code and 
the physical processes that we include. Then, we apply the 
code to perpendicular and oblique shocks with parameters 
similar to those of DRD including a single-sized grain fluid 
(Sect. [3]) or multiple grain fluids (Sect. [4}. Finally, we discuss 
the results and give our conclusions in Sect. [5] 



2 THE MODEL 
2.1 Numerical code 

As the ionisation fraction within molecular clouds is low, the 
plasma needs to be treated as a multicomponent fluid con- 
sisting of neutrals, ions, electrons and N dust grain species. 
The governing equations for the neutral fluid are given by 

dp n , d , s \ - 

dt 
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where r / / is the resistivity along the magnetic field, m the 
Hall resistivity and tad the ambipolar res istivity. Th e ex- 
pressions of the resistivities are given by e.g. lFallel (|2003l ) and 
depend strongly on the Hall parameter, /3j = on\B\/K n ip n , 
which is the ratio between the gyrofrequency of the charged 
fluid and collision frequency of the charged fluid with the 
neutrals. Substitution of the above expression into the 
Maxwell-Faraday equation then gives the magnetic induc- 
tion equation 
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where M = v n x B is the hyperbolic magnetic flux and R 
is the resistance matrix which depends on the Hall resistiv- 
ity, the ambipolar resistivity, and the resist ivity a long t he 
magnetic field. (For an expression for R, see iFallel (2003).) 

Equations [T| - [Tj are solved using the numerical scheme 
described bv lFalld (2003). This scheme uses a second-order 
hydrodynamic Godunov solver for the neutral fluid equa- 
tions (Eqs. IH3J1 . The charged fluid densities are calculated 
using an explicit upwind approximation to the mass conser- 
vation equation (Eq. [4| , while the velocities and tempera- 
tures are calculated from the reduced momentum and energy 
equations (Eq. [5] and [6}. The magnetic field is advanced ex- 
plicitly, even though this implies a rest riction on t he stable 
time step at high numerical resolution (Falle 2003). 
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where e n = p n /{"/ — 1) + is the internal energy per 

unit mass, s jn the mass transfer rate between the charged 
fluid j and the neutral fluid, Hj the energy sinks or sources 
for the different fluids, B the magnetic field, E the electric 
field and J the current given by J = OjPjVj. 

In the limit of small mass densities for the charged flu- 
ids, their inertia can be neglected. Also, the charged fluids 
are in thermal balance as their heat capacities are small. For 
the charged fluid j, the equations then reduce to 



The numerical code described in the previous section is quite 
general. In order to investigate the effects of dust grains 
in fast-mode shocks it is necessary to specify the source 
terms, i.e. the mass, momentum and energy transfer coef- 
ficients and the energy sources and sinks, to describe the 
relevant physics. The chemical network that we ado pted 
is the one used by IPilipp. Hartquist. fc Havnesl (|l990l ) and 
iPilipp fc Hartauistl i| 19941 ) . 
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where a.j is the charge to mass ratio, Kj„ the collision rate 
coefficient of the charged fluid j with the neutrals and Gj n 
the energy transfer rate per unit volume between the charged 
fluid and th e neutral fluid . The expressions for Kj„ and Gj n 
are given bv lDraindl| 19861 ). We use the reduced energy equa- 
tion to calculate the temperatures of the ion and electron 
fluids. For the grain fluids, a hydrodynamic approximation 
with zero pressure is appropriate as the thermal velocity 
dispersion of the dust grains is small compared to the drift 
velocity with the neutrals. 

The collisions of the charged particles with the neutral 
particles affect the evolution of the magnetic field. Using the 
reduced momentum equations for the charged fluids (Eq. [SJ) 
and the expressi on for the cur rent J, the electric field can 
be expressed as (|Cowlindll95r3 ) 



Mass transfer between the ion and electron fluids and 
the neutral fluid is due to cosmic ray ionisation at a rate 
of 10 -17 s _1 , electron recombination with Mg + , dissociative 
recombination with HCO + and recombination of ions on 
charged grains. Here we have a ssumed that the main gas 
phase ions are Mg+ and HCO+ (lOppenh eimer fc Dalgarnq 
Il974i ). The cross-section for these recombination reactions 
are taken from IPilipp. Hartquist. fc Havnesl (Il990l ). 

Mass is also transferred between the ion and electron 
fluids and the grain fluids due to ion neutralisation on grain 
surfaces. Ions and electrons stick to dust grains upon col- 
liding with them. Also, they pass on their charge with unit 
efficiency. The grain charge then critically depends on the 
electron-grain and ion-grain collision rates (which them- 
selves depends on the electron density and temperature, 
the ion density and temperature, and the ion-grain rela- 
tive streaming speed). The average grain charge is evalu- 
ated following lHavnes. Hartquist. fc Pilippl (|l987l ). Individ- 
ual grains can have a charge that differs from the aver- 
age value, but it is mostly within one unit of the average 
<|Elmegreenll 19791 ). 

Different radiative cooling processes for the neutral 
and electron fluids are included in the code. The elec- 
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tron gas is cooled by the excitation and subsequent de- 
excitation of molecular hydrogen (e.g. DRD). This process 
only becomes important when the electron fluid temper- 
ature exceeds » 1000 K. In the neutral gas, the domi- 
nant cooling processes are associated with O, CO, H2 and 
H2O. The radiative transitions for atomic oxygen, i.e. 01 
lines, are calculated using the procedure of DRD. The ra- 
diative losses from excited rotational-vibrational states of 
molecular hydrogen are a fit to the data presented by 
lHartauist. Dalgarno. fc Qppenheimei (Il980h. while t hose o f 
CO are calcula ted by using iHollenbach fc McKeei l| 19791 ). 
Finally, we use iNeufeld fc Melnickl (|l987f T to calculate the 
losses due to rotational transitions in H2O. The heating rate 
of the neutral gas by cosmic rays is from DRD. We have not 
included any radiative losses for the ion fluid . From Eq. [6] 
and using the expression for d i|Drainelll986h we find that 
the ion temperature, Ti, is given by 



Akb 



Vi) 
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where T n is the neutral temperature, m n the neutral particle 
mass and fcs the Boltzmann constant. 



2.3 Initial conditions 

The initial conditions of our simulatio ns are a d opted from 
DRD and similar studies such as IWardlei |l998) and 
IChapman fc Wardlei (|2006h . This allows a comparison of our 
results with those of these studies. 

The neutral fluid has a number density of nu — 
10 4 cm -3 or TiH = 10 6 cm -3 and is composed mainly of 
molecular hydrogen so that m n ~ 2mu- Small amounts 
of other neutral species such as O, CO, and H2O are also 
present. The initial abundances of O, CO, and H2O are taken 
to be 4.25xl0 -4 tih, 5 x 10~ 5 riH and respectively. How- 
ever, as the abundances of O and H2O evolve substantially as 
gas is heated, the abundances of the important coolants are 
recalculated at every time-step. The dominant ions within 
the gas are assumed, as mentioned above, to be HCO + and 
Mg + . The ion mass is then rtii ~ 30mH- Also, we assume 
that the ions are only singly ionised. 

The dust grains are all assumed to be spherical with a 
radius of a g = OAfim and a mass of 8.04 xlO" 13 g (where 
the density of the grains is taken 3 g cm -3 ). The grain to hy- 
drogen mass ratio in the unshocked region has the canonical 
value of 0.01. The densities of the ion and electron fluids and 
the average grain charge are calculated from ionisation equi- 
librium. For the temperatures that we adopt, i.e. T n = 8.47^ 
for n H = 10 4 cm" 3 and T„ = 26.7K for n H = 10 6 cm" 3 , the 
average grain charge is quite close to — e and the electron 
and ion number densities are roughly the same, i.e. n e ~ rii. 

In the interstellar medium the magnitude of the mag- 
netic field scales roughly as B ps l/LtG(riH/cm~ 3 )~ 1 ' /2 (see 
DRD) so that B = 10~ 4 G for n H = 10 4 cm" 3 and B = 10" 3 
G for «h = 10 6 cm -3 . The Alfven speed of the gas is then 
the same at both densities, i.e. va = 2.2 km s -1 . In the far 
upstream region the magnetic field lies in the x, i/-plane and 
makes an angle of 30° , 45° , 60° or 90° (perpendicular shock) 
with the z-axis. We assume that a shock is propagating in 
the positive x-direction with a speed of 25 km s _1 . Such a 
shock is a strong fast-mode shock with an Alfvenic Mach 
number of Ma ~ 11.5. 



Using the upstream fluid properties, the properties far 
downstream of the shock are calculated from the Rankinc- 
Hugionot relations for an isothermal, magnetised shock. For 
our initial shock structure we assume that the shock is a 
discontinuity at x = separating the upstream and down- 
stream properties. We follow the evolution of the shock to 
steady state within the shock frame. 



2.4 Computational details 

The size of the computational box is assumed to be a few 
time the shock width. By balancing the ion-neutral drag 
force with the Lorentz force, we ca n estimate the shock 
thickness (e.g. iDraine fc McKeell 19931 ) 
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where v s is the shock speed. However, the grain-neutral fric- 
tion is comparable to the ion-neutral friction for 

pg Ki 

or 
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where we substitut ed the expres sions for the collision rates 
with the neutrals l|Draine|[l986T ) and assumed the relative 
speed between the neutrals and the grains to be of the order 
v s /2. This constraint is fulfilled in the heated precursor for 
both densities and the shock thickness is determined by the 
grain-neutral drag. The expression for the shock width is 
similar to Eq. [9] and is given by 



(10) 



Values of y appropriate for the precu r sor m ust be used. 
iFlower. Pineau des Forets. fc Hartquistl (|l985l ) first pointed 
out that chemistry leads to a large drop in \ in the precur- 
sor of many shocks. The numerical domain is then chosen to 
be — 10L 9 ^ x ^ 10L g with a free-flow boundary condition 
at the downstream side and fixed inflow at the upstream 
boundary. We use a uniform grid with 400 cells. 

As our initial conditions are discontinuous, the inertia 
of the charged fluid, even the ions, cannot be neglected dur- 
ing the initial stages of the evolution. However, the inertial 
phase of the charged fluid is short com pared to the time 
to ap proach a steady configuration (e.g. iRoberge fc Ciolekl 
120071 ). Although the inertial phase does not affect the final 
steady shock structure here, this effect should not be ignored 
for some other time-dependent models. 



3 C-TYPE SHOCK MODELS 

In this section we present steady C-type shock structures 
for large single-sized grains. First we discuss perpendicular 
shocks and then shocks propagating at an oblique angle with 
the magnetic field. 
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Figure 1. The shock structure for a perpendicular shock propagating in a quiescent region with njj = 10 4 cm — 3 (left) and njj = 10 6 cm -3 
(right). From top to bottom we plot the velocity in the x-direction and z-direction for the neutrals (solid), ions/electrons (dashed) and 
grains (dotted), the tangential magnetic field and the temperatures of the ions (dashed), electrons (dotted) and the neutrals (solid) 
as well as the absolute value of the grain charge (dash-dotted). The velocities are normalised with respect to the shock speed and the 
magnetic field to the total upstream magnetic field. 



3.1 Perpendicular shocks 

Figure Q] shows the structures in multiple flow variables for 
perpendicular shocks at different densities. It is clear that 
both shock structures represent C-type shocks. As our mod- 
els are similar to DRD, we can compare our results with 
figs. 2 and 3 of DRD. Although the shapes of the shocks are 



qualitatively similar, our shock widths are about 4 times 
larger than in DRD. This is due to a lower fractional ionisa- 
tion in our models which is a consequence of including the 
chemistry governing the ionisation rather than assuming a 
constant ion flux as DRD did. The drag force of the charged 
particles on the neutrals is then lower as well and is bal- 
anced by the magnetic pressure over a larger length scale. A 
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Figure 2. The shock structure for an oblique shock propagating at an angle of 30° (left), 45° (middle) and 60° (right) with the magnetic 
field in a quiescent region with rin = 10 4 cm -3 . From top to bottom we plot the velocity in the ^-direction, y-direction and ^-direction for 
the neutrals (solid), ions/electrons (dashed) and grains (dotted), the tangential magnetic field in the y-direction (solid) and ^-direction 
(dashed), and the temperatures of the ions (dashed), electrons (dotted) and the neutrals (solid) as well as the absolute value of the 
grain charge (dash-dotted). The velocities arc normalised with respect to the shock speed and the magnetic field to the total upstream 
magnetic field. 
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consequence of the larger shock widths is that the gas has 
more time to radiate away its energy. The maximum tem- 
perature of the neutrals within the shock front is only 530 
K for n H = 10 4 cm" 3 and 905 K for n H = 10 6 cm" 3 , while 
DRD found 1223 K and 1334 K respectively. 

Figure [T] shows that the speed of the grains relative to 
the neutrals is less than the speed of the ions relative to 
the neutrals. The difference in relative speeds is due to the 
different response of the charged fluids to the Pedersen along 
E' (= E + v n x B) and Hall current along E'xB within 
the shock front. From Eq. [5] we find 
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1+1% 



E'xB 

B 2 



+ 



y+m 



(ii) 



where /3j is the Hall parameter of the charged fluid. For 
electrons and ions with \j3 e ^\ S> 1, the particles are tied to 
the magnetic field and 

V e .l - v„ E'xB 



c B 2 

At the other extreme, \/3j \ <C 1, the charged particles move 
with the neutrals. For the nu — 10 6 cm -3 models, the grains 
have a Hall parameter ranging from -0.1 (the upstream 
value) to a minimum of -1.5 within the shock front. For 
these values, the E' term has a non-negligible contribution 
to the drift velocity. For the a;-component of the velocity 
of the grains, this term only becomes important when the 
grains switch from moving with the neutrals to moving with 
the ions and electrons (see Fig. [lj . The ^-component of the 
E' term, however, has a larger contribution to the drift forc- 
ing the grains to move with the neutrals. The lower density 
model exhibits this effect to a lesser extent as the grain 
Hall parameter lies between -1 and -20. Then the E' term 
is less important and provides only a small contribution to 
the drift. 

The velocity difference between the ions and electrons 
with the neutrals determines the temperature of the charged 
particles through collisional heating. The maximum velocity 
difference in the models is « v s giving a maximum ion tem- 
perature of a few times 10 4 K. Electron cooling becomes im- 
portant at temperatures of about 1000 K reflected in a con- 
siderably lower electron temperature at roug hly 3xl0 3 K. 
The electron temperature within the shock front is impor- 
tant for the dynamics of the shock as the average grain 
charge is 



Z g e. 



4khT e a a 



(12) 



for \Z g \ » 1 and n g \Z g \ < n e fe.g. |Prainelll980l ). Otherwise, 
the average grain charge is roughly — 1. From Fig. [1] we find 
that the average grain charge follows this dependence on the 
electron temperature with the average grain charge reaching 
maximum values of about -250. 



3.2 Oblique shocks 

For a finite magnetic field component parallel to the shock 
propagation and a non-zero Hall conductivity, the Hall cur- 
rent along E' x B has a component parallel to the propa- 
gation direction of the shock. As the total current in that 
direction is zero if the shock is steady, the Hall current needs 



to be balanced by the Pedersen current (along E'). This Ped- 
ersen current gives rise to the ro tation of the magnetic fiel d 
around the shock propagation fPilip p~fc Hartauistl ll994T ). 
Significant rotation occurs when the Hall conduct ivity is 
larger than the Pedersen conductivity (IWardlelll998h . 

For the tih = 10 4 cm -3 models, the Hall conductiv- 
ity is always smaller than the Pedersen conductivity and 
there is no significant rotation of the magnetic field (see 
z-component of magnetic field in Fig. [2}. The rotation in- 
creases marginally as the angle, 9, between the direction 
of shock propagation and the magnetic field decreases. The 
other shock properties, also, do not change significantly with 
9. The temperatures of the different fluids and the average 
grain charge are similar to the ones found for the perpendic- 
ular shock. These results agree well wi th the ^-dependence 
of the steady shock structure found bv lChapman fc Wardlei 
(2006). The similarity of the shock structure for different val- 
ues of 9 results from the minor contribution of the grains to 
the drag force on the neutrals. The grains do not determine 
the dynamics and drift with the other charged species. 

Significant differences, however, are seen in the nn = 
10 6 cm -3 models. Here, the grains are the dominant con- 
tributors to the drag force on the neutrals. Furthermore, 
there is a significant charge separation giving rise to a large 
Hall conductivity. In Fig. [3] we see that there is a substan- 
tial rotation of the magnetic field. Also the magnetic field 
upstream is no longer a stable node, but a spiral node as 
the Hall conductivity cha nges wave propagation upstream 
|Wardlelll99Sl ; iFallel 120031 ). As for n H = 10 4 cm" 3 , the rota- 
tion decreases for larger 9, while the shock width increases. 

The rotation of the magnetic field has large conse- 
quences for the drift velocity of the charged particles. Fig- 
ure [3] shows that the velocity structure changes consider- 
ably with 9. Also, the grains move significantly differently 
than the ions and electrons, especially in the 9 = 30° model 
(see Fig. [3j . Along the propagation direction, electrons and 
ions are decelerated immediately. The grains, however, are 
initially accelerated before being decelerated. Furthermore, 
they lag the ions and electrons. To drift in the direction op- 
posite to the ions and electrons, the grains motion must be 
dominated by the E' term in Eq. [TT] and E' x must be in the 
same direction as the a;-component of E' x B. The large dif- 
ferences in drift velocity are expected as the rotation of the 
magnetic field is induced by a charge separation (and thus 
large Hall conductivity). For larger 9, the Hall conductivity 
and rotation decrease, so that the ions, electrons, and grains 
move more together. As expected, the velocity profile con- 
verges to the perpendicular shock structure (see Fig. [T] and 

El. 



4 MULTIPLE GRAIN SPECIES 

In Sect. [3] we have assumed a single-sized grain fluid. In 
molecular clouds , however, the dust has a grain-siz e distri- 
bution given by (M athis. Rumpl. fc Nordsieckll 19771 ) 



dn 
da 



(13) 



where dn is the number density of grains with radii between 
a and a + da. To take into account the effect of multiple 
dust species, we include an additional small-grain fluid with 
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Figure 3. Same as in Fig. [2] but with an upstream density of npj = 10 6 cm 



a„ — 0.04/im (i.e. a s = a g /W) and m 3 — 8.03 x 10 _16 g (or 
m 3 = m 9 /1000). The total mass density of the grain fluid, 
p s + p g , is assumed to be 1% of the neutral mass density (as 
in our previous models). With the adopted grain-size distri- 
bution most of the mass is in the small grain fluid. Further- 
more, the small grains dominate the grain-neutral friction 



and, thus, the shock dynamics as the collision coefficient be- 
tween grains and neutrals is proportional to Og j3 /m 3iS (see 
lDrainelTl986t ). From Fig. 2J we see that the shock width 
is indeed, as expected from Eq. If 01 an order of magnitude 
smaller than for the model with single-sized large grains. 
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Although the small grains have an average charge about 
an order of magnitude smaller than the large grains (see 
Eg. I12p . the small grains have a larger Hall parameter and 
are strongly coupled to the magnetic field. The velocity plots 
in Fig. [4] show that the small grains move with the electrons 
and ions, while the large grains move with velocities between 
those of the neutrals and the other charged particles. The 
grain charge density is dominated by the small grains and 
free electrons are depleted from the gas phase (i.e. rii 3> n e ). 
The charge separation between the charged particles is thus 
small and the Hall conductivity is too. The rotation of the 
magnetic field out of the plane is then negligible. Also, the 
upstream magnetic field is no longer a spiral node, but a 
stable node. 

The inclusion of a second grain species thus changes 
the shock structure considerably. Including additional grain 
species will further affect the shock structure as more dust 
grains are loading the magnetic field lines. Then the signal 
speed at which disturbunces propagate through the charged 
fluid decreases. For small PAH abundances the signal speed 
is close to 25 km s" 1 (|Flower fc Pineau des Foretsl I2003T ) 
leading to the possibility that the resulting shock is a J-type 
shock in which the grain inertia cannot be neglected. 
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Figure 4. Similar to the 45° model of Fig. \3\ but for two grain 
fluids. The small grains are indicated on the velocity components 
figures with a thin-dotted line, while their average charge is given 
by a thin dashed-dotted line. 



5 SUMMARY AND DISCUSSIONS 

In this paper we have presented ID time-dependent multi- 
fluid MHD simulations of perpendicular and oblique C-type 
shocks in dusty plasmas. We have included in our numerical 
code relevant mass transfer processes, such as electron re- 
combination with Mg + and dissociative recombination with 
HCO + , and radiative cooling by O, CO, H2 and H2O. Also, 
the mass and charge transfer from ions and electrons to the 
dust grains is taken into account to calculate the average 
grain charge. Our models are the first of oblique fast-mode 
shocks in dense molecular clouds in which a fluid treatment 
of grain dynamics has been combined with a self-consistent 
calculation of the thermal and ionisation balances including 
appropriate microphysics. We have performed simulations 
with single-sized and with two dust grains species in plas- 
mas of different density. 

Dust grains do not change the shock dynamics much at 
nu = 10 4 cm -3 and the shock structures for oblique and per- 
pendicular shocks are quite similar. Grains, however, domi- 
nate the shock dynamics at high densities (hh = 10 6 cm~ 3 ). 
For oblique shocks, the magnetic field in the shock front 
is not confined to the upstream magnetic field plane. The 
amount of rotation obtained depends on the Hall and Peder- 
sen conductivities. For single-sized grains the rotation is sig- 
nificant and increases with smaller 6 as the Hall conductivity 
in these models is comparable to the Pedersen conductivity. 
For some models, the magnetic field at the upstream side of 
the shock even spirals around the direction of shock propa- 
gation. When we include an additional small grain species, 
the rotation mostly disappears. The small grains which 
carry most of the charge move with the ions and electrons 
reducing the Hall conductivity. Our models corroborate 



Pilipp & Hartquist 1994; Wardle 1998; Cha 


oman & Wardle! 


2006; Guillet, Pineau des Forets. & Jones! 


20071), showing 



that our numerical method is efficient, rigorous, and robust. 
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In addition, we have succeeded in taking a significant step 
by developing the first oblique shock code capable of pro- 
ducing reliable results for comparisons with observations. 
Previous oblique shock models were restricted by less rigor- 
ous treatments of either the grain dynamics or the thermal 
and ionisation balances. 

However, our model does have some restrictions. Firstly, 
we determine the drift velocities of the grains solely by bal- 
ancing the grain-neutral friction with the Lorentz force and 
neglect their inertia. While small dust grains can be con- 
sidered inertialess, large grains have a drag length, i.e. the 
length scale for gas drag to decelerate the grains, that in 
some cases is comp arable to the length s cale for magnetic 
field compression (|Ciolek fe Robergd l2002). Therefore, large 
grain inertia cannot always be neglected. Secondly, for small 
grains, the average grain charge is close to the elementary 
charge. As there is some fluctuation around the average, a 
fraction of the small grains has a zero or even positive charge. 
The neutral grains decouple from the magnetic field, while 
the positive grains drift in a different direction to the neg- 
ative grains. Although this fluctuation around the average 
grain cha rge changes the shock structure, its effec t is only 
minimal dGuillet. Pineau des Forets. fe Jones! 120071 ). 

A time- dependent method is not only robust (even 
for non-equilibrium conditions) for finding steady 
state shock structures. It is also vital for modelling 
transient phenomena. In the vicinity of protostellar 
objects the abundance of SiO is significant l y en- 
hanced (e.g. Martin-Pintado. Bachiller. fe Fuentd 1 19921 ; 



Ijimenez-Serra et al.ll2004h . It is thought that the SiO emis- 
sion arises from the interaction of a p rotostellar jet with 
local dense clumps (|Lefloch et al.l fl998l ) . Grain sputtering 
within C-type shocks rel eases the SiO depleted on the grains 
back into the gas p hase (jCaselli. Hartauist. fe Havnes Il997l : 
ISchilke et ai1ll997h . Numerical models using steady C-type 
shocks indeed show good agreement with recent obser- 
vations of SiO line intensities in the L 1157 and L1448 
molecular outflows (Gu sdorf et al.l 2008a, b). However, it is 
doubtful that shocks interaction with dense clumps have 
reached steady state. In subsequent papers, we will study 
the interaction of C-type shocks with dense clumps and 
include grain sputtering to calculate the SiO emission. 
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